Cosmology with gravitational lensing 
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Abstract In these lectures I give an overview of gravitational lensing, concentrating 
on theoretical aspects, including derivations of some of the important results. Topics 
covered include the determination of surface mass densities of intervening lenses, 
as well as the statistical analysis of distortions of galaxy images by general inhomo- 
geneities (cosmic shear), both in 2D projection on the sky, and in 3D where source 
distance information is available. 3D mass reconstruction and the shear ratio test are 
also considered, and the sensitivity of observables to Dark Energy is used to show 
how its equation of state may be determined using weak lensing. Finally, the article 
considers the prospect of testing Einstein's General Relativity with weak lensing, 
exploiting the differences in growth rates of perturbations in different models. 



1 Introduction 

Gravitational lensing has emerged from being a curiosity of Einstein's General Rel- 
ativity to a powerful cosmological tool. The reasons are partly theoretical, partly 
technological. The traditional tool of observational cosmology, the galaxy redshift 
survey, has become so large that statistical errors are very small, and questions about 
the fundamental limitations of this technique were raised. Fundamentally, studies of 
galaxy surveys will be limited by uncertain knowledge of galaxy formation and 
evolution - we will probably never know with high accuracy where galaxies should 
exist, even statistically, in a density field. Since it is the density field which is most 
directly predictable from fundamental theories (with some caveats), gravitational 
lensing is attractive as it is a direct probe. Furthermore, statistical analysis shows 
that large weak lensing surveys covering large fractions of the sky to a median 
redshift of around unity are in principle extremely powerful, and can lead to er- 
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ror bars on cosmological parameters which are very small. In particular, the subtle 
effects of Dark Energy and even modifications to Einstein's General Relativity are 
potentially detectable with the sort of surveys which are being planned for the next 
decade. The limitations of lensing are likely to be in systematic errors arising from 
the ability to measure accurately shapes of galaxy images, physical effects aligning 
galaxies themselves, uncertain distances of sources, and uncertainties in the theoret- 
ical distribution of matter on small scales where baryon physics becomes important. 
Specially-designed instrumentation and telescopes with excellent image quality help 
the situation a lot, and understanding of the physical systematics is now quite good. 
These notes will not be concerned very much with the practical issues, but more 
with the theoretical aspects of how lensing works, how it can be used to determine 
surface densities of clusters of galaxies (testing the Dark Matter content), how it can 
be used statistically on large scales (testing the Dark Energy properties) and how 
the combination of geometrical measurements and the growth rate of perturbations 
can probe the gravity law (testing so-called Dark Gravity). It is not a comprehen- 
sive review, and in particular does not cover the observational developments, which 
have seen the typical size of lensing surveys increase from 1 (past) 100 (now) 
— >■ 10"* (near future) square degrees. For more comprehensive recent reviews, see 
e.g. ll62l[3Tl or the excellent SAAS-FEE lecture notes |65|. The structure of these 
notes is that section 2 covers basic lensing results, section 3 (Dark) Matter mapping, 
section 4 lensing on a cosmological scale, section 5 3D lensing and Dark Energy, 
and section 6 Dark Gravity. 



2 Basics of lensing 



We begin with some basic results on light bending, and derive results for the bending 
of light by an intervening lens with an arbitrary surface density. 



2.1 The bend angle 



The basic mechanism for gravitational lensing is that a point mass M will deflect a 
light beam through an angle 

4GM 

a{r) = — ^ Point Mass (1) 

where G is Newton's gravitation constant, c is the speed of light, and r is the distance 
of closest approach of the ray to the mass. This deflection angle is calculated using 
Einstein's General Theory of Relativity (GR), but is simply a factor two larger than 
Newtonian theory would predict if one treats the photon as a massive particle. The 
light path is curved, but if we are looking at a source which is far beyond the lens, the 
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light path can be approximated by two straight Hnes, with a bend angle a between 
them. This is the thin-lens approximation. 

The bend angle for a point mass equation ([TJ can be used to calculate the bend 
angle for an arbitrary mass distribution. For example, the bend angle for a thin lens 
whose mass distribution projected on the sky is circularly-symmetric depends on 
the projected mass enclosed within a circle centred on the middle of the lens, and 
extending out to the projected radius R of the ray; 

a(R) = — — ^ — ^ Circularly — symmetric mass (2) 
where M(< R) is the mass enclosed within a projected radius R. 



2.2 The lens equation 

The lens equation relates the true position on the sky of a source to the position of its 
image(s). We can get the basic idea of gravitational lensing by considering a point 
mass, or a circularly-symmetric mass distribution. Fig.[T]shows the geometry of the 
situation. The lens equation comes simply by noting that PSi = PS + SS\ (where P 
is the unmarked point directly above L). If we denote the observer-source distance 
by Ds, the lens-source distance by D^s, the position on the sky of the image by an 
angle 9, and the position on the sky of the unlensed source by /3, then, for small 
angles, this translates to 

Dse=Dsl5+DLsa (3) 

and we note that the bend angle is a function of the distance of the impact parameter 
of the ray R, a — a{R — DiQ), where Dl is the distance to the lens. For cosmolog- 
ical applications, the D, are angular diameter distances. 
For convenience, we define the scaled bend angle by 

Dls 

a=—a. (4) 
Ds 

This gives us the lens equation 

j3 = - a. (5) 

Note that for a given source position /3, this is an implicit equation for the im- 
age position(s) 9, and for general circularly-symmetric mass distributions, the lens 
equation cannot be solved analytically. There may indeed be more than one solu- 
tion, leading to multiple images of the same source. The lens equation is, however, 
an explicit equation for the source position given an image position. This is a use- 
ful feature which can be exploited for ray-tracing simulations of thin-lens systems, 
when one can ray-trace backwards from the observer to the sources. 
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Fig. 1 The geometry for a thin lens (adapted from J. Wambsganss). 
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2.2.1 Point mass lens: Multiple images and Einstein Rings 

The generic possibility of multiple images can be illustiated nicely by a point mass 
lens, where the implicit lens equation can be solved analytically for the image posi- 
tions. For a point mass M, the bend angle is a{R) = 4GM/{Rc^) = AGM/{Dl6c'^), 
so the lens equation is 

= (6) 



c20 DlDs 



which is a quadratic for Q : 



02-/30-02=0 (7) 
where we have defined the Einstein Angle 

There are evidently two images of the source for a point mass lens, one either side 
of the lens, at angles of 

0± = I ± + (9) 



The solutions are the intersections of the line Oi{d) and the straight line — j3. For 
< 0, the bend is in the opposite direction, so a < 0. We see that, with the point 
mass, there are inevitably two solutions for j3 7^ 0. 

There is a special case when the observer, lens and source are lined up (/3 = 0), 
when one singular solution is the straight-line path, and the other, at = is an 
Einstein Ring. Clearly, from the rotational symmetry of the arrangement, the image 
will appear at 9e in all directions round the lens. A near-perfect Einstein ring is 
shown in the colour plate section. 

Notice that if |j3| Be, then the main image is perturbed only slightly, at j3 + 
0J//3 +0{ji^^), and the second image is very close to the lens at 0_ ~ — 0j/j3. As 
we will see shortly, this second image is very faint. Be is a useful ballpark angle for 
determining whether the deflection of the source rays is significant or not. Typically, 
it is about an arcsecond for lensing by large galaxies at cosmological distances, and 
microarcseconds for lensing of stars in nearby galaxies by stars in the Milky Way. 



2.2.2 Magnification and Amplification 

As we have seen, simple lenses alter the positions of the image of the source, and 
may indeed produce multiple images. Another important, and detectable, effect is 
that the apparent size and the brightness of the source will change when its light un- 
dergoes a lensing event. The gravitational bending of light preserves surface bright- 
ness, so the change in apparent size is accompanied by a similar change in bright- 
ness. 
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For the circularly-symmetric lenses, the ratio of the solid angle of the source and 
that of the image gives the amplification of an infinitesimally small source: 

e de 

^ = m- '''' 

For a point lens, the amplifications of the two images are obtained straightforwardly 
by differentiating equation (j9]): 
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Note that the amplification can be negative. This corresponds to an image which is 
flipped with respect to the source. 

In some cases, such as microlensing, where the two images are unresolved, one 
can only measure the total amplification, 

Atotai = |A+ 1 + |A_ I = ^, ^ " (12) 

and we note that the difference of the amplifications of images by a point lens is 
unity. In the limit |j3 1 > 0£, A+ 1 + (Oe/P)'^ and A_ -{Oe/P)'^, so the inner 
image is extremely faint in this Umit. 



2.3 General thin lens mass distributions 



For a general distribution of mass, the bend angle is a 2D vector on the sky, a. The 
lens equation then generalises naturally to a vector equation 

I3 = e-a. (13) 

For a thin lens with surface density L{0'), the bend due at to a small element 
of mass dM — E{0') d^9' is evidently in the direction — & , and has a magnitude 
AGdM/{c^\e-e'\), so the total bend angle is 

r.(a\ ^GDeDes f ,2^, :E{0'){0-0') .... 

The bend angle can be expressed as the (2D) gradient of the (thin lens) lensing 
potential, a. = V0, where 

^{9) = '^^^^^ j d^0'Z{e')\n{\e~e'\) (tWnlens). (15) 
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This potential satisfies the 2D Poisson equation 

V20(0) =2k-(6») (16) 
where V is a 2D gradient, with respect to angle. K is the convergence k{0) = 

i:(0)/Zerit, With 

c^Ds 

^crit = ^ (17) 

AtzGDlDls 

being the critical surface density. The 2D Poisson equation is extremely useful, as 
it allows us to estimate the surface mass density of an intervening lens from lensing 
measurements, as we shall see later. 



2.3.1 Convergence, magnification and sliear for general tliin lenses 

For a general distribution of mass within a thin lens, the magnification and distortion 
of an infinitesimal source are given by the transformation matrix from the source po- 



Au^l^^S<j-^ij^ (18) 



sition /3 to the image position(s) 0. From the vector lens equation (13 1. The (inverse) 
amplification matrix is 

where = d^(j) /dOidOj. We see that A is symmetric, and it can be decomposed 
into an isotropic expansion term, and a shear A general amplification matrix also 
includes a rotation term (the final degree of freedom being the rotation angle), but 
we see that weak lensing doesn't introduce rotation of the image, and has only 3 
degrees of freedom, rather than the four possible in a 2 x 2 matrix. We decompose 
the amplification matrix as follows: 

^-Co^".)-(:^rf) 

where K is called the convergence and 

r=ri+'72 (20) 

is the complex shear. For weak lensing, both \k\ and \Yi\ are <C 1. A non-zero K 
represents an isotropic expansion or contraction of a source; 7i > represents an 
elongation of the image along the x-axis and contraction along y. 71 < stretches 
along y and contracts along x. 72 7^ represents stretching along x — ±.y directions. 
The effects are shown in Fig|2] 

Making the decomposition, we find that 
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=^(011+022) (21) 

71 = ^(011 -te) =Oi0 

72 = 012 = £>20- 

which defines the two operators D12, and it is straightforward to prove that 

D.Di = (V2)2, (22) 

where the summation is over i— 1,2. 

Note that K > corresponds to magnification of the image. Lensing preserves 
surface brightness, so this also amounts to amplification of the source flux. The 
magnification is 

^"d^" (l-)f)2-|7|2- ^^^^ 

We see that we may have infinite amplifications if fc > 1 . Such effects apply only for 
infinitesimal sources, and places in the source plane which lead to infinite magnifi- 
cations are referred to as caustics, and lead to highly distorted images along lines 
called critical lines in the lens plane. The giant arcs visible in images of some rich 
clusters lie on or close to critical lines. For cosmic shear, due to the general inho- 
mogeneities along the line-of-sight, K and |7| are typically 0.01, and the lensing is 
weak. 

It is worth noting that the amplification matrix may be written 

where g = 7/(1 — if) is called the reduced shear. Since the 1 — if multiplier affects 
only the overall size (and hence brightness) of the source, but not its shape, we see 
that shear measurements can determine only the reduced shear, and not the shear 
itself. For weak lensing, if <C 1, so the two coincide to linear order 

Note that a rotation of n/4 in the coordinate system changes a real shear into an 
imaginary shear - i.e. the complex shear rotates by n /2, twice the angle of rotation of 
the coordinate system. This behaviour is characteristic of a spin-weight 2 field, and 
is encountered also in microwave background polarisation and gravitational wave 
fields. 

The shear field is in principle observable, so let us see how it can be used to 
estimate the surface mass density of a lens. 



2.3.2 Estimating shear 

The estimation of shear from galaxy image data is a complex business, and 1 will 
give no more than a highly simplified sketch. For more details, see reviews such as 
ll62ll79ll3TI . One way to estimate shear is to measure the complex ellipticity e of a 
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Fig. 2 Ttie effect of shear dis- 
tortions on a circular source. 
In the notation of the current 
paper, e, and are the real 
and imaginary parts of the 
ellipticity (or shear). From 
(79J. 




galaxy, which can be defined in terms of moments ||46ll even if the galaxy image is 
not elliptical in shape. Fig. |2] shows how simple shapes map onto ellipticity. In the 
limit of weak distortions, the observed ellipticity is 

e = ^ (25) 

l + fe' 

where e" is the undistorted source ellipticity. The source ellipticity has a dispersion 
(|e'P) ~ 0.3, whereas yitself is usually much smaller, 0.01-0.1. If we average over 
some sources, 

{e) = Y (26) 

so we can use the average ellipticity as an estimator of shear 
Two points to note: 

1 . The estimator is noisy, dominated by the intrinsic dispersion of e, so 7 has a 
variance ^ /N if sources are averaged. 

2. In these notes the formulae will often refer to 7; it should be noted that in 
reality it will be e, a noisy estimate of 7, which is used in practice. 

The main practical difficulty of lensing experiments is that the atmosphere and 
telescope affect the shape of the images. These modifications to the shape may arise 
due to such things as the point spread function, or poor tracking of the telescope. The 
former needs to be treated with great care. Stars (whose images should be round) 
can be used to remove image distortions to very high accuracy, although a possibly 
fundamental limitation may arise because of the finite number of stars in an image. 
Interpolation of the anisotropy of the PSF needs to be done carefully, and examples 
of how this can be done in an optimal way are given in II8OI . Bayesian methods are 
beginning to be used, with lensfit 16111511 showing very promising results. 
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Theoretical work with numerical simulations indicates that in the absence of the 
effects of baryons, virialised Dark Matter haloes should follow a uniform 'NFW 
profile, p(r) = p,{rjr){l + r/r,)-^ 1631, if the Dark Matter is cold (CDM). Sim- 
ulations also predict how the physical size of the clusters should depend on mass, 
characterised by the concentration index q = ryir/vs, where r,,„ is defined as the 
radius within which the mean density is 200 times the background density. Roughly 
C.S °^ M^** ' . This can be tested, by measuring the shear signal and stacking the results 
from many haloes to increase signal-to-noise. 



3.1 2D mass surface density reconstruction 



We wish to take a map of estimated shear (ellipticities) and estimate the surface mass 
density of the intervening lens system. For simplicity here we assume the sources 
are at the same distance. The classic way to do this was given by fl3]| . 
In some respects it is easiest to work in Fourier space. Expanding 

Kk = J f-eK{9)exp{ik.9) (27) 

etc, then evidently 

K-k = -^^^0k (28) 

rik = -^{kj - kl)(j)k 
72k = -kik2^ 



where =^kj+kl ^ \k\^. 

We see that the following are estimators of fCt: 

k'^- k^ 
k\ — k2 2kik2 

In practice, we use the shapes of galaxies to estimate shear, so our estimates of 7k 
are dominated by noise in the form of scatter in the intrinsic shapes of galaxies. Thus 
the variance of each of these estimators is determined by the variances in |e,kP- The 
variance of the first estimator is therefore proportional to k^/{k\ — k2f', and the 
second to k'^ / [2kik2f'- The optimal estimator for )fk is given by the standard inverse 
variance weighting, giving 
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/t2 



rik- 



72k- 



(30) 



Since this is the sum of two terms, each of which is a muhiplication in k space, it 
represents a convolution in real space. (An exercise for the enthusiastic reader is to 
Fourier transform to find the real-space convolution). In fact it is easier to find the 
convolution by noting that 



(31) 



Poisson equation - it is given in equation 15 

1 



where the last line follows from equation |22| Now we know the solution to the 2D 
give 

I d^e'ri{o')in\e'-e\ 



In 



Upon differentiation with D,- and summation, we find 



2 f^2g/ [yi(^0co s(2v/) + 72(00 sin(2y)] 



(32) 



(33) 



where y/ is the angle between 6 and 6'. 

It is quite tempting to work with equation 33 in its discrete form, replacing the 
integral by a sum over shear estimates at the galaxy locations g, i.e. considering the 
following estimator for k: 



k{0) ^ 1 ^ [7i(eg)c os(2y) + 72(g,)sin(2v/)] 



(34) 



where n is the mean surface density of sources. This estimator would, however, be 
a mistake. It is unbiased, but it has the awkward and undesirable property of having 
infinite noise. This can be seen as follows. 



The fourier transform of equation 34 
1 



'm ''2 

k^ 



£7i(e,)exp(/k.0,) 



2kik2 

k^ 



^72Wexp(/k.0,) 



(35) 



If we assume (not quite correctly - see a later discussion on intrinsic alignments; but 
this will do here) that the galaxy shapes are uncorrected, then 



2n 



(36) 



where (|e^|) is the dispersion in galaxy elUpticities. Note that as usual, the estimate 
of the shear is dominated by the intrinsic ellipticity. 
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If we use Parseval's theorem, we see that the variance in the recovered conver- 



gence, being the integral of equation 36 over k, diverges. 

The solution is to smooth the field, which can be done at any stage - either 
smoothing the estimated shears (or averaging them in cells), or by working in 
Fourier space and introducing a filter, or even by smoothing the estimated noisy 
convergence field. The amount of smoothing required can prevent high-resolution 
maps of all but the richest clusters, but the advent of high surface densities with 
space telescopes has improved this situation considerably. 

Finally, we note that the estimate of the zero-wavenumber coefficient K"o is unde- 
termined by equation 35 We cannot determine the mean surface density this way. 
This feature is an example of the so-called 'mass-sheet degeneracy'. Other measure- 
ments (amplification/magnification, or by observing sources at different distances 
behind the lens) can alleviate this problem, or one can assume that if — > as one 
goes far from the lens centre. 



Fig. 3 2D reconstruction of 
surface matter density from 
the C0MB017 A901 study 
(contours; |21 1) superimposed 
on galaxy surface density. 




3.2 Testing the Navarro-Frenk-White profile of CDM 

Fig. |4] shows the average radial surface density profiles for clusters identified in the 
Sloan Digital Sky Survey (SDSS), grouped by number of cluster galaxies, and NFW 
fits superimposed |58 1. Fig|5]shows that the observed concentration indices are close 
to the theoretical predictions, but some tension exists. Broadly, weak lensing data 
on clusters therefore supports the CDM model. We will consider 3D mapping later, 
but unfortunately the limited accuracy of photo-zs (^ 0.03 ~ lOOh^ ^ Mpc typically) 
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Fig. 4 Excess surface density 
from stacked galaxy clusters 
from the SDSS survey, with 
best-fitting NFW profiles. 
^200 is a measure of the 
richness of the clusters. From 

m 

means the 3D mass map is smoothed heavily in the radial direction, and this limits 
the usefulness of 3D mapping for testing the NFW profile. 

A more radical test of theory has been performed with the Bullet Cluster ifTSll . 
actually a pair of clusters which have recently passed through each other There 
are two clear peaks in the surface density of galaxies, and X-ray emission from hot 
shocked gas in between. In the standard cosmological model, this makes perfect 
sense, as the galaxies in the clusters are essentially collisionless. If the (dominant) 
Dark Matter is also collisionless, then we would expect to see surface mass con- 
centrations at the locations of the optical galaxy clusters, and this is exactly what is 
observed. In MOND or TeVeS models without Dark Matter, one would expect the 
surface mass density to peak where the dominant baryon component is - the X-ray 
gas. This is not seen. A caveat is that it is not quite the surface density which is 
observed, rather the convergence, which is related to the distortion pattern of the 
galaxy images, and which is proportional to surface density in General Relativity 
(GR), but not in MOND/TeVeS. However, no satisfactory explanation of the bullet 
cluster has been demonstrated without Dark Matter. 



4 Cosmological lensing 

Gravitational lensing is strong if the distortions to the images are substantial, and 
weak if the distortions are small. Weak lensing of background images can be used 
to learn about the mass distribution in individual objects, such as galaxy clusters, 
but we concentrate now on weak lensing on a cosmological scale, which is usually 
analysed in a statistical way. Reviews of this area include 1621 [3T1 |6l l69l l79l . 




R [h-'Mpc] 
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Fig. 5 Concentration indices 
from SDSS clusters as a func- 
tion of mass, compared with 
simulation (dotted, for dif- 
ferent cosmologies). Dashed 
regions show range assuming 
a power-law cs — M relation. 
From [58j|. 

The basic effect of weak lensing is for the clumpy matter distribution to perturb 
slightly the trajectories of photon paths. By considering how nearby light paths are 
perturbed, one finds that the shapes of distant objects are changed slightly. Asso- 
ciated with the size change is a change in the brightness of the source. Size and 
magnitude changes can, in principle, be used to constrain the properties of the mat- 
ter distribution along the line-of-sight (and cosmological parameters as well), but 
it is the change in shape of background sources which has almost exclusively been 
used in cosmological weak lensing studies. The reason for this is simply that the 
signal-to-noise is better These notes will concentrate on shear (=shape changes), 
but the magnification and amplification of sources can also be used and will prob- 
ably be used in future when the surveys are larger The great promise of lensing is 
that it acts as a direct probe of the matter distribution (whether dark or not), and 
avoids the use of objects which are assumed to trace the mass distribution in some 
way, such as galaxies in large-scale structure studies. Theoretically, lensing is very 
appealing, as the physics is very simple, and very robust, direct connections can be 
made between weak lensing observables and the statistical properties of the matter 
distribution. These statistical properties are dependent on cosmological parameters 
in a known way, so weak lensing can be employed as a cosmological tool. The main 
uncertainties in lensing are observational - it is very challenging to make images of 
the necessary quality. In this section, we concentrate here on the weak effects on a 
cosmological scale of the non-uniform distribution of matter all the way between 
the source and observer, an effect often refeiTed to as cosmic shear. 
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4.1 Distortion of light bundles 



The distortion of a light bundle has to be treated with General Relativity, but if one 
is prepared to accept one modification to Newtonian physics, one can do without 
GR. 

In an expanding Universe, it is usual to define a comoving coordinate x, such 
that 'fundamental observers' retain the same coordinate. Fundamental observers are 
characterised by the property of seeing the Universe as isotropic; the Earth is not 
(quite) a fundamental observer, as from here the Cosmic Microwave Background 
looks slightly anisotropic. The equation of motion for the transverse coordinates 
(about some fiducial direction) of a photon in a flat Universe is 

d^Xi 2 d4> . , ^ 

We assume a flat Universe throughout, but the generalisation to non-flat Universes 
is straightforward (there is an extra term in the equation above, and some r symbols 
need to be changed to an angular diameter distance). 

The equation of motion can be derived using General Relativity (see Appendix 
for details). We will use jc, , / = 1 , 2 for coordinates transverse to the line-of-sight, and 
r to indicate the radial coordinate. Tj is the conformal time, related to the coordinate 
thydr\= cdt/R{t) and R{t) is the cosmic scale factor, equal to Rq at the present 
time, ^'(x,, r) is the peculiar gravitational potential, related to the matter overdensity 
field 5 = 8p/p by Poisson's equation 

(38) 

where Hq is the Hubble constant, Q,„ is the present matter density parameter, and 
a{t) =R{t)/Ro = (1 where z is redshift. 

The equation of motion is derived in the Appendix from General Relativity, in a 
(nearly flat) metric given in the Newtonian gauge by 

ds^ = (1 +2<I>/c^)c^dt^ - {1 -2(I>/c^)R^{t){dr^ + r^de^ + r^sin^ ed(p^). (39) 

From a Newtonian point-of-view, equation[37]is understandable if we note that time 
is replaced by rj (which arises because we are using comoving coordinates), and 
there is a factor 2 which does not appear in Newtonian physics. This same factor of 
two gives rise to the famous result that in GR the angle of light bending round the 
Sun is double that of Newtonian theory. 

The coordinates x, are related to the (small) angles of the photon to the fiducial 
direction = (0^, 0,.) by x, = 
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4.2 Lensing potential 



The solution to ( 37 1 is obtained by first noting that the zero-order ray has ds^ = 
c/r = —drj, where we take the negative root because the light ray is incoming. 
Integrating twice, and reversing the order of integration gives 

Xi = rOi~^ / dr'^{r-r'). (40) 
Jo dXj 

We now perform a Taylor expansion of (9^>/(9x; , and find the deviation of two 
nearby light rays is 

Axi^rAe,-^Ae, [dr'r\r-r')-^^ (41) 







which we may write as 

Axi = rAej{5ij-^ij) (42) 
where 5/; is the Kronecker delta (/ = 1 , 2) and we define 

,„„,^r.,(:zP^. ,43) 

Jo rr dBidBj 

The integral is understood to be along a radial line (i.e. r || r'); this is the Born 
approximation, which is a very good approximation for weak lensing |[8] |66l Ell • In 
reality the light path is not quite radial. 

It is convenient to introduce the (cosmological) lensing potential, which controls 
the distortion of the ray bundle: 

2 /■'■ ir — r') 

(j){r) = ^ dr , — ^(r') (cosmological; flat universe) (44) 

c^ Jo rr 

Note that 0,j(r) — d^<p{r) / dOidOj. So, remarkably, we can describe the distortion 
of an image as it passes through a clumpy universe in a rather simple way. The shear 
and convergence of a source image is obtained from the potential in the same way as 
the thin lens case (equation 22 1, although the relationship between the convergence 
and the foreground density is more comphcated, as we see next. 



4.2.1 Relationship to matter density field 

The gravitational potential ^> is related to the matter overdensity field 5 = 5p/p by 
Poisson's equation ( [38] l. The convergence is then 

.(r)^5|^r*'^^|^ (45) 

Ic^ Jo r a[r') 
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Note that there is an extra term d'^<t>/dr'^ in which integrates to zero to the 
order to which we are working. 



4.2.2 Averaging over a distribution of sources 

If we consider the distortion averaged over a distribution of sources with a radial 
distribution p{r) (normalised such that / dr p{r) = 1), the average distortion is again 
obtained by reversing the order of integration: 

where 

g{r)^ Tdr'pir')'^. (47) 

J r f 

In order to estimate p{r), surveys began to estimate distances to source galaxies 
using photometric redshifts. This has opened up the prospect of a full 3D analysis 
of the shear field, which we will discuss briefly later in this article. 



4.2.3 Convergence power spectrum and Shear correlation function 

The average shear is zero, so the most common statistics to use for cosmology are 
two-point statistics, quadratic in the shear These may be in 'configuration' ('real') 
space, or in transform space (using Fourier coefficients or similar). I will focus on 
two quadratic measures, the convergence power spectrum and the shear-shear cor- 
relation function. 

To find the expectation value of a quadratic quantity, it is convenient to make 
use of the matter density power spectrum, P{k), defined by the following relation 
between the overdensity Fourier coefficients: 

(5k5k*,) = (27r)-^5^(k-k')P(fc), (48) 

where 5^ is the Dirac delta function. P{k) is evolving, so we write it as P{k;r) in 
future, where r and t are related through the lookback time. (This r-dependence may 



look strange; there is a subtlety: (48 i holds if the field is homogeneous and isotropic, 
which the field on the past light cone is not, since it evolves. In the radial integrals, 
one has to consider the homogeneous field at the same cosmic time as the time of 
emission of the source). The trick is to get the desired quadratic quantity into a form 
which includes P{k; r). 

For the convergence power spectrum, we first transform the convergence in a 2D 
Fourier transform on the sky, where £ is a 2D dimensionless wavenumber: 
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Ki = / d-'9K{e)e 



-ite 



drr 



a{r) 



d^e5{re,r)e 



(49) 
(50) 



where A = 3Hq£2,„/2c^. We expand the overdensity field in a Fourier transform, 



d{re,r) 



{2n) 



(51) 



and substitute into (50 1. We form the quantity {KfK"^,), which, by analogy with (48 1, 
is related to the (2D) convergence power spectrum by 

{K,K}) = {2nf8''{l-i')PM- (52) 

Straightforwardly, 



drG{r) / dr'G{r') / d^Od^O' 



, d'k d^k' 



{2n)^ (27r)3 

(5k5^/)exp(/A:||r-/;t||V)exp(/kx.6»-/k'L.0')exp(-/^.6> + /£'.6>') 



(53) 



where G{r) = rg{r)/a{r). Using (48 1 we remove the k' integration, introducing the 
power spectrum =P{^kJ^ + |kxp). For small-angle surveys, most of the signal 

comes from short wavelengths, and the k^^ is negligible, so P{k:) ~ P(|kx|). The 
only term remaining is the exponential, which integrates to (2;r)5^(r — r'). The 
integrals over 9 and 0' give (2;r)^5^(£ — rkx) and {27i)^d^{i' — rkj_') respectively, 
so the whole lot simplifies to give the convergence power spectrum as 



3//q i2„, 
2c2 



dr 



a{r) 



P{£/r;r). 



(54) 



An exercise for the reader is to show that the power spectrum for the complex shear 
field 7 is the same: Py = P^- The shear correlation function, for points separated by 
an angle Q is 



(rf)6 



d^i 
m 
dl 



jte 



(55) 



"fd(p 



2% 



fMe)jQ{£e) 



where we have used polar coordinates, with (p the angle between £ and 6, and we 
have exploited the isotropy (P^ depends only on the modulus of £). Jo is a Bessel 
function. 
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Other quadratic quantities (examples are shear variances on different scales, 
Aperture Mass (squared)) can be written similarly as integrals over the power spec- 
trum, with different kernel functions. 



4.3 Matter power spectrum 

As we have seen, the two-point statistics of the shear and convergence fields de- 
pend on the power spectrum of the matter, P{k;t). The power spectrum grows in 
a simple way when the perturbations in the overdensity are small, |5| ^ 1, when 
the power spectrum grows in amplitude whilst keeping the same shape as a func- 
tion of k. However, gravitational lensing can still be weak, even if the overdensity 
field is nonlinear. Poisson's equation still holds provided we are in the weak-field 
limit as far as General Relativity is concerned, and this essentially always holds for 
cases of practical interest. In order to get as much statistical power out of lensing, 
one must probe the nonlinear regime, so it is necessary for parameter estimation to 
know how the power spectrum grows. Through the extensive use of numerical sim- 
ulations, the growth of dark matter clustering is well understood down to quite small 
scales, where uncertainties in modelling, or uncertain physics, such as the influence 
of baryons on the dark matter |81|, make the predictions unreliable. Accurate fits 
for the nonlinear power spectrum have been found |71 1 up to A: > lOhMpc^ ' , which 
is far beyond the linear/nonlinear transition k ^ O.lhMpc^^ . Fig. |6] shows fits for 
a number of CDM models. For precision use, one must make sure that the statis- 



Fig. 6 The nonlinear power 
spectrum from numerical 
simulations, along with fitting 
functions (from L71i . 
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tics do not have substantial contributions from the high-A: end where the nonlinear 
power spectrum is uncertain. This can be explored by looking at the kernel functions 



4.4 Intrinsic alignments 

The main signature of weak lensing is a small alignment of the images, at the level 
of a correlation of ellipticities of ^ 10^*. One might be concerned that physical pro- 
cesses might also induce an alignment of the galaxies themselves. In the traditional 
lensing observations, the distances of individual galaxies are ignored, and one sim- 
ply uses the alignment on the sky of galaxies, and one might hope that the galaxies 
will typically be at such large separations along the line of sight that any physical 
interactions would be rare and can be ignored. However, the lensing signal is very 
small, so the assumption that intrinsic alignment effects are sufficiently small needs 
to be tested. This was first done in a series of papers by a number of groups in 
2000-1 e.g.|l23l [TSl [161 [131 , and the answer is that the effects may not be negligi- 
ble. The contamination by intrinsic alignments is highly depth-dependent. This is 
easy to see, since at fixed angular separation, galaxies in a shallow survey will be 
physically closer together in space, and hence more likely to experience tidal inter- 
actions which might align the galaxies. In addition to this, the shallower the survey, 
the smaller the lensing signal. In a pioneering study, the alignments of nearby galax- 
ies in the SuperCOSMOS survey were investigated [11 1. This survey is so shallow 
(median redshift ^0.1) that the expected lensing signal is tiny. A non-zero align- 
ment was found, which agrees with at least some of the theoretical estimates of the 
effect. The main exception is the numerical study of Jing |44|, which predicts a 
contamination so high that it could dominate even deep surveys. For deep surveys, 
the effect, which is sometimes called the TI correlation', is expected to be rather 
small, but if one wants to use weak lensing as a probe of subtle effects such as the 
effects of altering the equation of state of dark energy, then one has to do something. 
There are essentially two options - either one tries to calculate the intrinsic align- 
ment signal and subtract it, or one tries to remove it altogether. The former approach 
is not practical, as, although there is some agreement as to the general level of the 
contamination, the details are not accurately enough known. The latter approach is 
becoming possible, as lensing surveys are now obtaining estimates of the distance 
to each galaxy, via photometric redshifts (spectroscopic redshifts are difficult to ob- 
tain, because one needs a rather deep sample, with median redshift at least 0.6 or 
so, and large numbers, to reduce shot noise due to the random orientations of ellip- 
ticities). With photometric redshifts, one can remove physically close galaxies from 
the pair statistics (such as the shear correlation function) 1261 l48l . Thus one removes 
a systematic error in favour of a slightly increased statistical error The analysis in 
ll27l is the only study which has explicitly removed close pairs. 

Another effect which is more subtle is the correlation between the orientation or 
a foreground galaxy and the orientation of the lensed image of a background galaxy. 
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The latter is affected gravitationally by the tidal field in the vicinity of the former 
(as well as the tidal fields all the way along the line-of-sight), and if the orientation 
of the foreground galaxy is affected by the local tidal field, as it surely must at 
some level, then there can be a contamination of the cosmological lensing signal by 
what is sometimes referred to as the 'GI correlation' . This was first pointed out in 
[30], and it seems to be likely to be a significant effect |29|. It is less easy to deal 
with than the 11 correlation, but modelling and nulling methods exist, at the price of 
diminished signal-to-noise lfT0ll40l . 



Fig. 7 Example patterns from 
E-mode and B-mode tields 
(from 1 79 1). Weak lensing 
only produces E-modes at 
any significant level, so the 
presence of B-modes can 
indicate systematic errors. 
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4.5 E/B decomposition 

Weak gravitational lensing does not produce the full range of locally linear distor- 
tions possible. These are characterised by translation, rotation, dilation and shear, 
with six free parameters. Translation is not readily observable, but weak lensing 
is specified by three parameters rather than the four remaining degrees of freedom 
permitted by local affine transformations. This restriction is manifested in a number 
of ways: for example, the transformation of angles involves a 2 x 2 matrix which 



is symmetric, so is not completely general, see equation ( 18 i. Alternatively, a gen- 
eral spin-weight 2 field can be written in terms of second derivatives of a complex 
potential, whereas the lensing potential is real. There are many other consistency 
relations which have to hold if lensing is responsible for the observed shear field. In 
practice the observed ellipticity field may not satisfy the expected relations, if it is 
contaminated by distortions not associated with weak lensing. The most obvious of 
these is optical distortions of the telescope system, but could also involve physical 
effects such as intrinsic alignment of galaxy ellipticities, which we will consider 
later. 
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A convenient way to characterise the distortions is via E/B decomposition, where 
the shear field is described in terms of an 'E-mode', which is allowed by weak 
lensing, and a 'B-mode', which is not. These terms are borrowed from similar de- 
compositions in polarisation fields. In fact weak lensing can generate B-modes, but 
they are expected to be very small 1671 . so the existence of a significant B-mode 
in the observed shear pattern is indicative of some non-lensing contamination. The 
easiest way to introduce a B-mode mathematically is to make the lensing potential 
complex: 

+ (56) 

There are various ways to determine whether a B-mode is present. A neat way is to 
generalise a common statistic called the aperture mass to a complex M ~ Map + /Mx, 
where the real part picks up the E modes, and the imaginary part the B modes. 
Alternatively, the ^± can be used llT7ll68l : 

PK±{i)--7t r dee[Ja{ie)^+{e)±Mie)^-{e)] (57) 

Jo 

where the ± power spectra refer to E and B mode powers. In principle this requires 
the correlation functions to be known over all scales from to °°. Variants of this 
(see e.g. f\n\) allow the E/B-mode correlation functions to be written in terms of 
integrals of ^± over a finite range: 

^E{e) = \[^-{e) + ^i{d)] (58) 

where 

^;(0) = ^+(0)+4y^ —Ue')-i29^l ^Ue')- (59) 

This avoids the need to know the correlation functions on large scales, but needs the 
observed correlation functions to be extrapolated to small scales; this was one of the 
approaches taken in the analysis of the CFHTLS data ll32ll . Difficulties with estimat- 
ing the correlation functions on small scales have led others to prefer to extrapolate 
to large scales, such as in the analysis of the GEMS [28 1 and William Herschel data 
ll59l . Note that without full sky coverage, the decomposition into E and B modes is 
ambiguous, although for scales much smaller than the survey it is not an issue. 



4.6 Results 

The first results from cosmic shear were published in 2000 ||3] ITT] |47] [82ll . so as 
an observational science, cosmological weak lensing is very young. To date, the 
surveys have been able to show clear detections of the effect, and reasonably ac- 
curate determination of some cosmological parameters, usually the amplitude of 
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the dark matter perturbations (measured by the rms fractional fluctuations in an 
Sh^^Mpc sphere and denoted cig), and the matter density parameter Q,„. Current 
surveys cannot Hft a near-degeneracy between these two, and usuafly a combina- 
tion (typically a^Qj^^) is quoted. This makes sense - it is difficult, but certainly not 
impossible, to distinguish between a highly-clumped low-density universe, and a 
modestly-clumped high-density universe. There is no question that the surveys do 
not yet have the size, nor the careful control of systematics required to compete 
with the microwave background and other techniques used for cosmological param- 
eter estimation. However, this situation is changing fast, particularly with the CFHT 
Legacy Survey, which is now complete, and Pan-STARRS 1, which is underway. 
Future surveys such as proposed by Euclid, JDEM and LSST should be much more 
powerful. Early results from CFHTLS are shown in Fig.|8]and|9] 



CQ 



Fig. 8 E- and B-modes from 
an early analysis of CFHTLS 
data (2). Top points are the 
E-modes and bottom the 
B-modes. 
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Fig. 9 Cosmological parame- 
ters from an early analysis of 
CFHTLS data fcTJ 
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5 Lensing in 3D 

Knowing the redshift distribution of sources is vital to interpret 2D lensing statistics, 
and to help with this, recent lensing surveys have obtained multi-colour data to 
estimate distances using photo-zs. With photo-zs for sources, it makes sense to use 
the information individually, and this opens up the possibility of 3D lensing. 



5.1 3D potential and mass density reconstruction 

As we have already seen, it is possible to reconstruct the surface density of a lens 
system by analysing the shear pattern of galaxies in the background. An interesting 
question is then whether the reconstruction can be done in three dimensions, when 
distance information is available for the sources. It is probably self-evident that 
mass distributions can be constrained by the shear pattern, but the more interesting 
possibility is that one may be able to determine the 3D mass density in an essentially 
non-parametric way from the shear data. 

The idea |73| is that the shear pattern is derivable from the lensing potential 0(r), 
which is dependent on the gravitational potential ^> (r) through the integral equation 



<.(r): 



dr' 



0(r') 



(60) 



where the integral is understood to be along a radial path (the Born approximation). 



and a flat Universe is assumed in equation (60 1. The gravitational potential is related 



to the density field via Poisson's equation (38 i. There are two problems to solve 
here; one is to construct ^ from the lensing data, the second is to invert equation 



(60 1. The second problem is straightforward: the solution is 



0(r) = 



2J~r 



r^^^Hr) 



(61) 



From this and Poisson's equation V^^> 
the mass overdensity field 



{3/2)HQ£2,„5/a{t), we can reconstruct 



5(r) 



dr 



(62) 



The construction of <j) is more tricky, as it is not directly observable, but must be 
estimated from the shear field. This reconstruction of the lensing potential suffers 
from a similar ambiguity to the mass-sheet degeneracy for simple lenses. To see 
how, we first note that the complex shear field 7 is the second derivative of the 
lensing potential: 
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r(r) 




0(r). 



(63) 



As a consequence, since the lensing potential is real, its estimate is ambiguous up to 
the addition of any field /(r) for which 



= 0. 



Since must be real, the general solution to this is 

/(r) = +G(r)0, +//(r)0,, +P(r)(02 + 02) 



(64) 



(65) 



where F, G, H and P are arbitrary functions of r = |r|. Assuming these functions 
vary smoothly with r, only the last of these survives at a significant level to the mass 
density, and corresponds to a sheet of overdensity 



5 = 



4fl(f)c2 d 



(•2 



(66) 



There are a couple of ways to deal with this problem. For a reasonably large survey, 
one can assume that the potential and its derivatives are zero on average, at each r, 
or that the overdensity has average value zero. For further details, see |4|. Note that 
the relationship between the overdensity field and the lensing potential is a linear 
one, so if one chooses a discrete binning of the quantities, one can use standard 
linear algebra methods to attempt an inversion, subject to some constraints such as 
minimising the expected reconstruction errors. With prior knowledge of the signal 
properties, this is the Wiener filter See [35| for further details of this approach. 

This method was first applied to COMBO-17 data \1A\, and recendy to COSMOS 
HST data iH - see FigfTUl 



5.2 Tomography 

In the case where one has distance information for individual sources, it makes 
sense to employ the information for statistical studies. A natural course of action 
is to divide the survey into slices at different distances, and perform a study of the 
shear pattern on each slice. In order to use the information effectively, it is necessary 
to look at cross-correlations of the shear fields in the slices, as well as correlations 
within each slice 1.33,1 . This procedure is usually referred to as tomography, although 
the term does not seem entirely appropriate. 

We start by considering the average shear in a shell, which is characterised by a 
probabihty distribution for the source redshifts z = z{r), p{z)- The shear field is the 
second edth derivative of the lensing potential llT2ll 
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Fig. 10 3D reconstmction 
of matter density from the 
COSMOS ACS data L6OI. 




(67) 



where the derivatives are in the angular direction, and the last equality holds in the 
flat-sky limit. If we average the shear in a shell, giving equal weight to each galaxy, 
then the average shear can be written in terms of an effective lensing potential 



MO) 



dzpiz)<j){r) 



(68) 



where the integral is at fixed 0, and p{z) is zero outside the slice (we ignore errors in 
distance estimates such as photometric redshifts; these could be incorporated with 
a suitable modification to p{z))- In terms of the gravitational potential, the effective 
lensing potential is 

0eff(0) = 4 rdr^{r)g{r) (69) 
Jo 

where reversal of the order of integration gives the lensing efficiency to be 

'1 1 



dz p(z') 
z{r] '^^ '\r r 



(70) 



z' = and we assume flat space. If we perform a spherical harmonic transform 
of the effective potentials for slices / and /', then the cross power spectrum can be 
related to the power spectrum of the gravitational potential P<p (k) via a version of 
Limber's equation: 

(71) 



where 



dr 



P<p{e/r,r) 



(72) 
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is the cross power spectrum of the lensing potentials. The last argument in f$ allows 
for evolution of the power spectrum with time, or equivalently distance. The power 
spectra of the convergence and shear are related to cffj by |I34T| 

Hij ^ H,7 ('^) 

yy _ 1 (^ + 2)! 

^ti] 4(^_2)!H<;- 

The sensitivity of the cross power spectra to cosmological parameters is through 
various effects, as in 2D lensing: the shape of the linear gravitational potential power 
spectrum is dependent on some parameters, as is its nonlinear evolution; in addition 
the z(r) relation probes cosmology, via 

riz)=cr^. (74) 

The reader is referred to standard cosmological texts for more details of the depen- 
dence of the distance-redshift relation on cosmological parameters. 




Hu f33\ illustrates the power and limitation of tomography, with two shells (Fig. 
[TT| ). As expected, the deeper shell (2) has a larger lensing power spectrum than the 
nearby shell (1), but it is no surprise to find that the power spectra from shells are 
correlated, since the light from both passes through some common material. Thus 
one does gain from tomography, but, depending on what one wants to measure, the 
gains may or may not be very much. For example, tomography adds rather little to 
the accuracy of the amplitude of the power spectrum, but far more to studies of dark 
energy properties. One also needs to worry about systematic effects, as leakage of 
galaxies from one shell to another, through noisy or biased photometric redshifts, 
can degrade the accuracy of parameter estimation ll36l l57l [U |50l |52l . 
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5.3 The Shear Ratio test 



The shear contributed by the general large-scale structure is typically about 1 %, but 
the shear behind a cluster of galaxies can far exceed this. As always, the shear of a 
background source is dependent on its redshift, and on cosmology, but also on the 
mass distribution in the cluster. This can be difficult to model, so it is attractive to 
consider methods which are decoupled from the details of the mass distribution of 
the cluster Various methods have been proposed (e.g. BTH gl). The method currently 
receiving the most attention is simply to take ratios of average tangential shear in 
different redshift slices for sources behind the cluster. 

The amplitude of the induced tangential shear is dependent on the source redshift 
Z, and on cosmology via the angular diameter distance-redshift relation 5jt[r(z)] by 

m 

where „ is the shear which a galaxy at infinite distance would experience, and 
which characterises the strength of the distortions induced by the cluster, at redshift 
Zi- Evidently, we can neatly eliminate the cluster details by taking ratios of tangential 
shears, for pairs of shells in source redshift: 



„ _ Yr, ^ Sk[r{zi)]Sk[rizi)-rizi)] 
''-yr, Sk[rizt)]S,[riZj)-rizi)y 



In reality, the light from the more distant shell passes through an extra pathlength 
of clumpy matter, so suffers an additional source of shear This can be treated as a 
noise term ll75l . This approach is attractive in that it probes cosmology through the 
distance-redshift relation alone, being (at least to good approximation) independent 
of the growth rate of the fluctuations. Its dependence on cosmological parameters 
is therefore rather simpler, as many parameters (such as the amplitude of matter 
fluctuations) do not affect the ratio except through minor side-effects. More signif- 
icantly, it can be used in conjunction with lensing methods which probe both the 
distance-redshift relation and the growth-rate of structure. Such a dual approach can 
in principle distinguish between quintessence-type dark energy models and mod- 
ifications of Einstein gravity. This possibility arises because the effect on global 
properties (e.g. z(r)) is different from the effect on perturbed quantities (e.g. the 
growth rate of the power spectrum) in the two cases. The method has a signal-to- 
noise which is limited by the finite number of clusters which are massive enough to 
have measurable tangential shear. In an all-sky survey, the bulk of the signal would 
come from the 10^ — 10^ clusters above a mass limit of IO'^Mq. 
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5.4 Full 3D analysis of the shear field 

An alternative approach to take is to recognise that, with photometric redshift esti- 
mates for individual sources, the data one is working with is a very noisy 3D shear 
field, which is sampled at a number of discrete locations, and for whom the lo- 
cations are somewhat imprecisely known. It makes some sense, therefore, to deal 
with the data one has, and to compare the statistics of the discrete 3D field with 
theoretical predictions. This was the approach of 1221 [T2l l24l l49ll . It should yield 
smaller statistical errors than tomography, as it avoids the binning process which 
loses information. 

In common with many other methods, one has to make a decision whether to 
analyse the data in configuration space or in the spectral domain. The former, usually 
studied via correlation functions, is advantageous for complex survey geometries, 
where the convolution with a complex window function implicit in spectral methods 
is avoided. However, the more readily computed correlation properties of a spectral 
analysis are a definite advantage for Bayesian parameter estimation, and we follow 
that approach here. 

The natural expansion of a 3D scalar field {r,9,(j)) which is derived from a poten- 
tial is in terms of products of spherical harmonics and spherical Bessel functions, 
jl(kr)Y^{0). Such products, characterised by 3 spectral parameters {k,£,m), are 
eigenfunctions of the Laplace operator, thus making it very easy to relate the expan- 
sion coefficients of the density field to that of the potential (essentially via —k^ from 
the operator). Similarly, the 3D expansion of the lensing potential. 




where the prefactor and the factor of k are introduced for convenience. The expan- 
sion of the complex shear field is most naturally made in terms of spin-weight 2 
spherical harmonics 2^7' and spherical Bessel functions, since 7= \ dd<^, and 

r(r) = V2^Y. I dkYimkjeXkr) zFf (0). (78) 

The choice of the expansion becomes clear when we see that the coefficients of the 
shear field are related very simply to those of the lensing potential; 

= (79) 

The relation of the (p(„ (k) coefficients to the expansion of the density field is readily 
computed, but more complicated as the lensing potential is a weighted integral of 
the gravitational potential. The details will not be given here, but relevant effects 
such as photometric redshift errors, nonlinear evolution of the power spectrum, and 
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the discreteness of the sampling are easily included. The reader is referred to the 
original papers for details. 

In this way the correlation properties of the Ye,„{k) coefficients can be related to 
an integral over the power spectrum, involving the z{r) relation, so cosmological 
parameters can be estimated via standard Bayesian methods from the coefficients. 
Clearly, this method probes the dark energy effect on both the growth rate and the 
z{r) relation. 



5.5 Dark Energy with 3D lensing methods 

In this section we summarise some of the forecasts for cosmological parameter esti- 
mation from 3D weak lensing. We concentrate on the statistical errors which should 
be achievable with the shear ratio test and with the 3D power spectrum techniques. 
Tomography should be similar to the latter We show results from 3D weak lens- 
ing alone, as well as in combination with other experiments. These include CMB, 
supernova and baryon oscillation studies. The methods generally differ in the pa- 
rameters which they constrain well, but also in terms of the degeneracies inherent 
in the techniques. Using more than one technique can be very effective at lifting 
the degeneracies, and very accurate determinations of cosmological parameters, in 
particular dark energy properties, may be achievable with 3D cosmic shear surveys 
covering thousands of square degrees of sky to median source redshifts of order 
unity. 



Fig. 5.5 and 13 show the accuracy which might be achieved with a number of 
surveys designed to measure cosmological parameters. We concentrate here on the 
capabilities of each method, and the methods in combination, to constrain the dark 
energy equation of state, and its evolution, parametrised by lfT4l 

w(a) = ^ = wo + w„(l -fl) (80) 

where the behaviour as a function of scale factor a is, in the absence of a compelling 
theory, assumed to have this simple form, w = —1 would arise if the dark energy 
behaviour was actually a cosmological constant. 

The assumed experiments are: a 5-band 3D weak lensing survey, analysed either 
with the shear ratio test, or with the spectral method, covering 10,000 square degrees 
to a median redshift of 0.7, similar to the capabilities of a groundbased 4m-class 
survey with a several square degree field; the Planck CMB experiment (14-month 
mission); a spectroscopic survey to measure baryon oscillations (BAO) in the galaxy 
matter power spectrum, assuming constant bias, and covering 2000 square degrees 
to a median depth of unity, and a smaller z = 3 survey of 300 square degrees, similar 
to WFMOS capabilities on Subaru; a survey of 2000 Type la supernovae to z = 1.5, 
similar to SNAP'S design capabilities. 

We see that the experiments in combination are much more powerful than indi- 
vidually, as some of the degeneracies are lifted. Note that the combined experiments 
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Fig. 12 The accuracy ex- 
pected from the combination 
of experiments dedicated 
to studying Dark Energy 
properties. Marginal Iff, 2- 
parameter regions are shown 
for the experiments individ- 
ually and in combination. 
The supernova study fills the 
plot, the thin diagonal band is 
Planck, the near-vertical band 
is BAO, and the ellipse is the 
3D tensing power spectrum 
method. The small ellipse is 
from combined. From 1241 . 

appear to have rather smaller error bars than is suggested by the single-experiment 
constraints. This is because the combined ellipse is the projection of the product of 
several multi-dimensional likelihood surfaces, which intersect in a small volume. 
(The projection of the intersection of two surfaces is not the same as the intersection 
of the projection of two surfaces). The figures show that errors of a few percent on 
Wo are potentially achievable, or, with this parametrisation, an error of w at a 'pivot' 
redshift of z ~ 0.4 of under 0.02. This error is essentially the minor axis of the er- 
ror ellipses. These graphs include statistical errors only; systematic errors, from, for 
example bias in the photo-z distribution, are expected to degrade errors by ~ a/2 
1501. 



6 Dark Gravity 

In addition to the possibility that Dark Energy or the cosmological constant drives 
acceleration, there is an even more radical solution. As a cosmological constant, Ein- 
stein's term represents a modification of the gravity law, so it is interesting to con- 
sider whether the acceleration may be telling us about a failure of GR. Although no 
compelling theory currently exists, suggestions include modifications arising from 
extra dimensions, as might be expected from string-theory braneworld models. In- 
terestingly, there are potentially measurable effects of such exotic gravity models 
which weak lensing can probe, and finding evidence for extra dimensions would of 
course signal a radical departure from our conventional view of the Universe. 
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Fig. 13 As in Fig. |53] but 
with the shear ratio test as the 
lensing experiment. Super- 
novae fill the plot, Planck is 
the thin diagonal band, BAO 
the near-vertical band, and the 
shear ratio is the remaining 
45 degree band. The combi- 
nation of all experiments is in 
the centre. From 1751 . 




Weak lensing is useful as it probes not just the distance-redshift relation, but 
also the growth rate of perturbations (see, for example, equations [70| and [72]i. This 
is important, because in principle measures which use only the r{z) relation suffer 
from a degeneracy between a different gravity law and the equation of state of the 
contents of the Universe. To see this, first note that a modified gravity law will lead 
to some sort of Hubble relation i/(a). If we combine the Friedmann equation 

H^{a) + ^^ (81) 



and 



we find that 



3 



^ ipa^) ^ ~pa^ ^ -w{a)na^ (82) 
da ' 



w(fl) = — 



3 d\na 



Q,n{a) 



1 



(83) 



So we see that any modified gravity law can be mimicked, as far as the distance- 
redshift relation is concerned, by GR with an appropriate equation of state. 

To analyse other gravity laws, we consider scalar perturbations in the conformal 
Newtonian gauge (flat for simplicity), ds^ = a^{ri) [(1 +2^)dr\^ — (1 — 20)(ix^], 
where i// is the potential fluctuation, and the curvature perturbation, T] being the 
conformal time. Information on the gravity law is manifested in these two potentials. 
For example in GR and in the absence of anisotropic stresses (a good approxima- 
tion for epochs when photon and neutrino streaming are unimportant) ^ = y/. More 
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generally, the Poisson law may be modified, and the laws for y/ and may differ. 
This difference can be characterised ||T9l by the slip, CJ. This may be scale- and 
time-dependent: \if{k,a) — [I +n}{k,a)] ^{k,a), and the modified Poisson equation 
may be characterised by Q, an effective change in G |2|: 

- k^(^ = AnGa^p,AQ{k, a). (84) 



Different observables are sensitive to \j/ and (p in different ways [42]. For example, 
the Integrated Sachs-Wolfe effect depends on \j/ + ([>, but the effect is confined to 
large scales and cosmic variance precludes accurate use for testing modified gravity. 
Peculiar velocities are sourced by Lensing is sensitive to yf+cj), and this is the 
most promising route for next-generation surveys to probe beyond-Einstein gravity. 
The Poisson-like equation for + is 

-k\xif + ^)=2E^^8,„ (85) 
Za 



where Z = 2(1 +CT/2). For GR, Z = 1, CJ = 0. The DGP braneworld model 11201 has 
= 1, so mass perturbations deflect light in the same way as GR, but the growth rate 
of the fluctuations differs. Thus we have a number of possible observational tests, 
including probing the expansion history, the growth rate of fluctuations and the mass 
density-light bending relation. Future WL surveys can put precise constraints on E 
lEl, and on d (see Fig. [TO) |19J. 



By probing the growth rate as well as the expansion history, weak lensing can lift 
a degeneracy which exists in methods which consider the distance-redshift relation 
alone, since the expansion history in a modified gravity model can always be mim- 
icked by GR and Dark Energy with a suitable w(a). In general however the growth 
history of cosmological structures will be different in the two cases (e.g. Il53ll37]| . 
but see L54j ). 



6.1 Growth rate 

Whilst not the most general, the growth index 7 ll55l (not to be confused with the 
shear) is a convenient minimal extension of GR. The growth rate of perturbations in 
the matter density p,„, 5„, = 8p„,/p,„, is parametrised as a function of scale factor 
a{t) by 

^ ^ g{a) = exp I r ^ [n,My- 1] I , (86) 
a [Jq fl' '- -'J 

In the standard GR cosmological model, 7~ 0.55, whereas in modified gravity theo- 
ries it deviates from this value. E.g. the flat DGP braneworld model |20| has 7 ~ 0.68 
on scales much smaller than those where cosmological acceleration is apparent ||56l . 

Measurements of the growth factor can in principle be used to determine the 
growth index 7, and it is interesting to know if it is of any practical use. In contrast to 
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0.35 



Fig. 14 The projected 
marginal 68% and 95% likeli- 
hood contours for the slip, W, 
assuming B7 = CTo(l +z)^^, 
for WMAP 5-year data (blue), 
adding current weak lensing 
and ISW data (red). Yellow is 
mock Planck CMB data, and 
green adds weak lensing from 
a 20,000 square degree survey 

EH. 
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Fig. 15 Expected Bayesian 
evidence B vs. deviation of 
the growth index from GR, for 
a future WL survey + Planck 
1251 . If modified gravity is the 
true model, GR will still be 
favoured by the data to the left 
of the cusp. The Jeffreys scale 
of evidence L43J is labeled. 
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parameter estimation, this is an issue of model selection - is the gravity model GR, or 
is there evidence for beyond-Einstein gravity? This question may be answered with 
the Bayesian evidence, B |70|, which is the ratio of probabilities of two or more 



models, given some data. Following 112511 . Fig. 15 shows how the Bayesian evidence 
for GR changes with increasing true deviation of 7 from its GR value for a combina- 
tion of a future WL survey and Planck. From the WL data alone, one should be able 
to distinguish GR decisively from the flat DGP model at InZ? ~ 11.8, or, in the fre- 
quentist view, 5.4c7 |25 1. The combination of WL+P/ancfe+BAO +SN should be able 
to distinguish 5/ = 0.041 at 3.41 sigma. This data combination should be able to de- 
cisively distinguish a Dark Energy GR model from a DGP modified-gravity model 
with expected evidence ratio InZ? ~ 50. An alternative is to parametrise geometry 
and growth with two separate effective values of w, and look for (in)consistency 
Il39l|72l|76l[83|. 

One caveat on all of these conclusions is that WL requires knowledge of the non- 
linear regime of galaxy clustering, and this is reasonably well-understood for GR, 
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but for other models, further theoretical work is required. This has already started 
||54l . The case for a large, space-based 3D weak lensing survey is strengthened, as it 
offers the possibility of conclusively distinguishing Dark Energy from at least some 
modified gravity models. 



Survey 


V 


|lnfi| 




DES+P/afic/t+BAO+SN 


3.5 


1.28 


substantial 


DES+Planck 


2.2 


0.56 


inconclusive 


DES 


0.7 


0.54 


inconclusive 


PSl+P/ancjt+BAO+SN 


2.9 


3.78 


strong 


PSl+Planck 


2.6 


2.04 


substantial 


PSl 


1.0 


0.62 


inconclusive 


WLNc+Planck+BAO+SN 


10.6 


63.0 


decisive 


WLf/c+Planck 


10.2 


52.2 


decisive 


WLng 


5.4 


11.8 


decisive 



Table 1 The evidence ratio for the three weak lensing experiments considered with and without 
Planck, supernova and BAO priors. WL/^a is a next-generation space-based imaging survey such 
as proposed for DUNE or SNAP. z,„ is the median redshift, no the number of sources per square 
arcminute, and cr. is the assumed photometric redshift error. For completeness, we also list the fre- 
quentist signiticance VCT with which GR would be expected to be ruled out, if the DGP braneworld 
were the correct model. 



7 The Future 

The main promise of weak lensing in the future will come from larger surveys with 
optics designed for excellent image quality. Currently, the CFHTLS is the state-of- 
the-art, covering ^ 170 square degrees to a median redshift in excess of one. In the 
near future Pan-STARRS, VST and DES promise very small PSF distortions and 
large areal coverage, and in the far future LSST on the ground, and satellites such as 
Euclid or IDEM may deliver extremely potent lensing surveys. In parallel with these 
developments, the acquisition of photometric redshifts for the sources has opened 
up the exciting possibility of analysing weak lensing surveys in 3D. Each source 
represents a noisy estimate of the shear field at a location in 3D space, and this 
extra information turns out to be extremely valuable, increasing substantially the 
statistical power of lensing surveys. In particular, it can lift the degeneracy between 
(78 and Q„,, measure directly the growth of Dark Matter clustering |5 | and, more 
excitingly still, it represents a powerful method to measure the equation of state of 
Dark Energy 1221 1411 - surely one of the most important remaining questions 
in cosmology. In addition, photometric redshifts allow the possibility of direct 3D 
Dark Matter mapping Il73l l4l l74l . thus addressing another of the unsolved problems. 
Finally, there is the interesting prospect that we may be able to find evidence for 
extra dimensions in the Universe, from braneworld models, or other cosmological 
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models, from the effect they have on the gravity law and the growth rate of perturba- 
tions. Most excitingly, it seems that this may be within reach for ambitious lensing 
surveys planned for the next decade. 
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8 Appendix: the propagation of light through a 
weakly-perturbed universe 

8.1 The geodesic equation 

The geodesic equation governs the worldline (A = 0, 1,2,3) of a particle, and is 
readily found in textbooks on General Relativity (e.g. [38 1). It is 

d^x^ dx^" dx'' _^ ^g^^ 

dp^ dp dp 

where p is an affine parameter, and F is the affine connection, which can be 
written in terms of the metric tensor g^y as 

^^^--2^ X^ + ^-^l ^^^^ 
For weak fields, the interval is given by 

ds^^ (^l + ^^cW-(^l-^^R^{t)[dr'- + Sl{r)dl5^] (89) 

where is the peculiar gravitational potential, R{t) is the scale factor of the Uni- 
verse, and r,0,(p are the usual comoving spherical coordinates. The angle c/j3 = 

d9^ + sin^ dd(p^. Sk^r) depends on the geometry of the Universe, being given by 

sinhr, if < 0; 
sinr, if A: > 0. 

The curvature A; = — 1 , 0, 1 corresponds to open, flat and closed Universes respec- 
tively. We are interested in the distortion of a small light bundle, so we can con- 
centrate on a small patch of sky. If we choose the polar axis of the coordinate sys- 
tem to be along the centre of the light bundle, we can define angles 9^ = 9 cos (p 
and 9y = 9sin(p. For convenience we also use the conformal time, defined by 
dr\ — cdt/R{t), in place of the usual time coordinate. With these definitions the 
interval is more simply written as 

ds^^R^t) { (l + ^) dri'-[l - ^) [dr^ + Slir)id9^ + d9^)] | . 



The metric tensor for weakly-perturbed flat Friedmann-Robertson-Walker metric is 
then 
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R\t) 



/l+20/c2 \ 

-(l-20/c2) 

-r2(l-20/c2) 

V -r^{\-2^/c^) ) 



We are interested in how the angles of the ray, (0^, 0>) change as the photon moves 
along its path, responding to the varying gravitational potential. The unperturbed, 
radial, path is set by = ds^ ~ dr\^ — dr^, i.e. For a radial incoming ray. 



dr 
dt] 



= -1. 



With ^'^^ defined as the inverse of (so gfj.vg'^" = 5^ the affine connections 
are readily computed. 

The parametrized equation for Tj is required only to zero-order in 0, and reduces 

to 

where a dot here denotes d/dp. By choosing the unit of p appropriately, we find 

d^ _ 1 
Jp^W' 

We can also relate the radial coordinate to the conformal time, again to zero-order 
in 4>: The first-order equations governing 0, and 0,, are obtained from the geodesic 
equation, or by the variational methods (see for example, d'Inverno (1992), section 
7.6 OSl) 



^jcM dp ydxl" J 



= 0, 



where = [ds/dpf. Withx'' = 0, 



j2 d0 J 



50/ " 



02)- 



dp 



-IR^? ( 1 



20 



Br 



0. 



With the zero-order solutions for drj/dp and dr/drj, to first order this reduces sim- 
ply to 

t/20_^ 2de^ _ 2 d<t> 

dr]^ r drj c^r^ dO^ 

It is convenient to write this as an equation for the comoving displacement of the 
ray from a fiducial direction, 

Xi = rdi. i — 1,2 
and the equation for 0^ and a similar one for Qy simplify to 
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(fx 




where V here is a comoving transverse gradient operator {dx,dy). 

We see that the propagation equation for the displacement looks similar to what 
one would guess from a Newtonian point-of-view; the presence of 77 (instead of f) 
in the acceleration term on the left is a result of the expansion of the Universe and 
the choice of comoving coordinates. The right-hand side looks like the gradient of 
the potential, but is larger than the naive gradient by a factor of two. This is the same 
factor of two which leads to the classic result of General Relativity, famously tested 
by Eddington's 1919 solar eclipse observations, that the angle of light bending by 
the Sun is double what Newtonian theory predicted. 



